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Abstract 

We have studied the reflection delay time distribution from a one-dimensional 
electronically random system of finite length. We show that the average delay time 
is a non-monotonic function of the strength of the disorder and exhibits three qual- 
itatively different regimes. In the weak disorder limit the scattering is dominated 
by the ordinary Born scattering. For the intermediate strengths of disorder a well 
defined regime of multiple scattering emerges followed by the localization domain. 
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1 INTRODUCTION 



Wigner delay time, the energy derivative of scat- 
tering phase shift, is considered to be related to 
the time spent by the particle in the region of 
interaction [1,2]. Since its inception in 1955 the 
Wigner delay time has been a quantity of profound 
interest from fundamental as well as technological 
point of view [3-6] . The delay time statistics is in- 
timately connected with the dynamic admittance 
of microstructures [3]. Although the phase delay 
time alone is not sufficient to calculate capacitance 
of extended mesoscopic systems, a direct relation 
between the capacitance and the phase time does 
exist for systems which can be treated essentially 
as zero dimensional [6] (eg. a cavity or a quantum 
dot). The wave packet incident on the surface of a 
sample is backscattered after some delay. In ran- 
dom systems this results in a low temperature 1/ / 



type noise for the fluctuating surface currents due 
to non-cancellation of the instantaneous currents 
at the surface [7]. For a three dimensional sys- 
tem of arbitrary shape with an arbitrary number 
of incoming channels the density of states is pro- 
portional to the sum of dwell times for all incom- 
ing channels [8] . In case of a one dimensional dis- 
ordered sample connected to an infinite lead the 
dwell time is the same as the Wigner delay time. 
This delay time is directly related to the density 
of states. 



In the single channel case the distribution of the 
delay time for a disordered semi-infinite sam- 
ple was obtained earlier by using the invariant 
imbedding approach [7]. The stationary distribu- 
tion P s (t) for the dimensionless delay time r is 
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given by 



2 NUMERICAL TECHNIQUE 



= (e x */ 2 -l)(l + r 2 ) 



where A is proportional to the disorder induced 
localization length and the most probable value 
of r occurs at T mp = A/2. The long time tail of 
the above distribution scales as 1/r 2 . The average 
value of r is logarithmically divergent indicating 
the possibility of the particle traversing the infi- 
nite sample before being totally reflected, due to 
the resonances [9,10]. If the disordered region is 
semi-infinite or if it is terminated at one end by 
an infinite wall, the reflection coefficient will be 
unity, and the complex reflection amplitude will 
have the form R = e ie{E) . 

A recent study based on analytical work found the 
delay time distribution in the one-channel case to 
be universal (especially the long time tail is inde- 
pendent of the nature of disorder) [11] within ran- 
dom phase approximation (RPA). Soon after the 
universality of tail distribution beyond RPA was 
established numerically [12]. The universality of 
stationary distribution of delay time was further 
supported by another analytical work where the 
random potential is of a different kind [13]. So far 
in this field the localization has been studied ex- 
tensively. Recently Texier et. al. [9] have studied 
the statistics of delay time distribution for ballis- 
tic case for which sample size is much smaller than 
the localization length. In the present paper we 
analyse the statistical properties of delay time and 
study the emergence of three qualitatively differ- 
ent regimes of transport or scattering in a finite 
one-dimensional random system. These different 
domains are dominated by ordinary Born scatter- 
ing, multiple scattering and localization respec- 
tively. The crossover from one regime to another 
as the disorder strength is varied is studied. 



The following Hamiltonian, written in a tight- 
binding, one-band model is used to describe the 
motion of a quasiparticle on a random lattice: 

H = ]T e n \n){n\ + F(|n)(n + l| + |n) (n - 1|), (2) 

where |n), e n and V denote the non-degenerate 
Wannier orbital at site n, the site energy at the 
site n and the hopping matrix element connecting 
nearest neighbors separated by a unit lattice spac- 
ing respectively. We consider site diagonal disor- 
der where the site energies e n are assumed un- 
corrected random variables distributed uniformly 
over -W/2 to W/2 (P(e n ) = l/W). The L site 
(n = 1 to L) disordered ID sample is connected 
at one end to a perfect infinite lattice having all 
site energies zero. The well-known transfer-matrix 
method [14,15] is used to calculate the reflection 
amplitude r(E) = \r\e~ te( - E ^ and its phase 6(E) 
at two values of incident energy E = Eq±5E. The 
delay time is then calculated using the definition 
r = UdO/dE. 

Throughout our following discussion we consider 
the delay time r in a dimensionless form by multi- 
plying them with V and we set fi = 1. For all the 
results discussed below we have used Eq = and 
5E = 10~ 8 . We have confirmed the stability of 
our results with respect to choice of SE within the 
range 10~ 6 < 5E < 10~ 9 . For the calculation of 
various distributions and averages we have used 
atleast 10 6 realizations of the disordered sample. 



3 Results and discussion 

In Fig. 1 we show the plot of average value of delay 
time (r) versus the disorder strength for a sample 
of length L = 50. For very small values of W such 
that localization length £ > L (£ = 96 /W 2 [15]), 
we see that the (r) decreases with increasing W. 
This is the ballistic regime. Since the disorder is 
very weak i.e., £ 3> L Born scattering dominates. 
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Due to this the waves undergo mostly single scat- 
tering events. The (r) continues to fall upto a 
certain strength of disorder depending upon the 
length of the sample. At this strength of disor- 
der £ ~ 20L. We have verified this for samples of 
different lengths. Increasing W beyond this value 
we observe that (r) starts increasing with disor- 
der, i.e., the disorder becomes strong enough such 
that multiple scattering events start dominating. 
This leads to particles (or waves) spending more 
time in the sample before getting reflected. Con- 
sequently, the reflection delay time increases with 
increasing disorder. This trend continues until the 
disorder becomes strong enough such that £ ~ L. 
On increasing W beyond this, a crossover to lo- 
calization occurs. In this regime the most prob- 
able distance traversed by the particles is twice 
the localization length. The localization length de- 
creases with disorder strength as 96 /W 2 . There- 
fore, the most probable value of the reflection 
delay time, the time required to travel to-and- 
fro distance of £, also decreases and so does (r). 
However, from the inset shown in the figure we 
see that the root mean squared relative variance 

a rms = \J (t 2 ) — (t) 2 / (r) in this regime is much 
larger than unity. This indicates that the fluctua- 
tions in r in this regime are very large and hence 
r is a non-self-averaging quantity. Thus, from the 
behaviour of average reflection delay time (r) we 
are able to identify three qualitatively different 
regimes of wave propagation in one-dimensional 
disordered medium. 

In Fig. 2 is shown the distribution P{r) of reflec- 
tion delay time r from a sample of fixed length 
L = 50 for three values of disorder strength W in 
the ballistic regime. For W = the distribution is 
a delta peak at a value of r equal to the round trip 
delay time tq. For the weak disorder case the most 
probable value of r is the round trip delay time To- 
From the figure we observe that the distributions 
are asymmetric about the most probable value 
To- The distribution has higher weight for values 
t < tq and lower weight for values of r > To- 
The asymmetry is confirmed by looking at the 
third moment of r which is Xs = ((t — (t)) 3 ) 
shown in the inset of Fig. 2. The variation of Xs 



with W sheds some light on the behaviour of (r) 
with W in the ballistic regime. In this regime, as 
mentioned earlier, disorder is weak and much of 
the behaviour is dominated by single scattering 
events which occur anywhere in the sample before 
reaching endpoint. As the disorder strength is in- 
creased the probability of scattering at individual 
sites also increases. Thus, the probability of par- 
ticle encountering a scattering near the surface 
itself increases with increasing disorder strength. 
This enhances the weight for small r in -P(r). Al- 
though, the multiple scattering events contribute 
to the distribution for r > To, the change in their 
contribution with increasing disorder is not much 
compared to that of single scattering events. 

The distributions of delay time in the multiple 
scattering regime and the localized regime are 
shown in Fig. 3. For the distribution in the mul- 
tiple scattering regime we notice that the most 
probable delay time T mp is much lower than 
the round trip delay time To- The distribution 
is broader and has a tail. The appearance of a 
broad distribution and a long time tail indicates 
that in this regime the multiple scattering events 
are indeed dominating. In fact, it is due to the 
appearance of a tail that the (r) increases with 
increasing disorder strength W in spite of the 
reduction in T mp . For W > \/2 (£ < L) we 
crossover to the localized regime. The correspond- 
ing distributions (for W = 1.4 and W = 2.4) are 
shown in the inset of Fig. 3. The T mp in this case 
is proportional to the typical length traversed by 
the particles that is the localization length £ [7]. 
Since £ decreases with increasing W so does T mp . 
The distribution shows a long-time-tail due to 
the presence of resonances. For these resonant 
realizations the particle explores the whole of 
the sample before being totally reflected. Also, 
for this case one can notice that the distribution 
becomes narrow and weight at the tails becomes 
smaller with increasing disorder in contrast with 
the behaviour observed for the tails in the multi- 
ple scattering domain. 

In the ballistic regime the transport is dominated 
by few scattering events and since the disorder is 
weak the phase changes associated with these are 
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small. Consequently, the distribution of phase of 
the reflected wave shows (Fig. 4) a doubly peaked 
structure in this regime. The peaks occur at and 
2tt. In the multiple scattering regime the multi- 
ple scattering events provide ample opportunity 
for the phase to get randomized. This is reflected 
in the fact that the distribution shown in Fig. 4 
for W = 0.7 (multiple scattering) is fairly uni- 
form over the range to 2tt. For W = 8.0 we are 
in the strong localization regime and once again 
the double peak structure emerges in the phase 
distribution. In this case, though the double peak 
structure arises because £ <C L i.e., the scattering 
is so strong that most of the particles are reflected 
back from near the surface thus getting scattered 
only from a few sites. In the limit of W — > oo the 
peaks shift towards tt/2 and 3tt/2. 

In Fig. 5 we plot the average value of time delay 
(r) as a function of length L of the sample. As 
predicted by Texier et. al. [9] we see that the de- 
pendence is linear on L. In the inset of the figure 
we plot, on log-log scale, the second cumulant of 
r (t 2 ) — (t) 2 as a function of the length L of the 
sample. We find a L 3 - 56 dependence of the second 
cumulant of r on the length of the sample in con- 
trast to the L 3 behaviour predicted by Texier et. 
al. [9]. The difference, probably, can be attributed 
to the approximations involved in the analytical 
calculations or the constraint of averaging over fi- 
nite number of realizations in our numerical cal- 
culations. 



4 Conclusion 



In summary, our study based on the statistics 
of Wigner delay time distribution, shows that 
three qualitatively different regimes of trans- 
port or scattering exist in disordered finite one- 
dimensional systems. The systems exhibit a suc- 
cessive crossover from ballistic regime to multiple 
scattering regime to localization regime as the 
strength of disorder is varied. 
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Fig. 1. Plot of (r) versus disorder strength W for 
different lengths of the disordered segment. The 
disorder averaging is done over 10 7 samples. 



Fig. 3. Distribution of r from a sample of length 
L = 50 for different values of disorder strength 
W. For L = 50, W = 0.5 and W = 0.7 lie in the 
multiple scattering regime. Inset shows the distri- 
bution of r for W = 1.4 and W = 2.4 belonging 
to the localization regime. 
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Fig. 2. Distribution of r from a sample of length 
50 for three different values of disorder strength 
W for the ballistic regime ^ «L, The inset shows 
the behaviour of third moment of r with W. 




Fig. 4. Distribution of phase of reflected 
wave from the semi-infinite disordered sample 
of length L = 50 for different values of disor- 
der strength W in the different regimes namely, 
ballistic(IF = 0.1), multiple scattering (W = 0.7) 
and strong localization (W = 8.0). 
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Fig. 5. Plot of Average delay time (r) versus 
length of the sample L for disorder strength 
W = 0.2. In the inset is the plot of second cumu- 
lant of r ((r 2 ) — (r) 2 ) versus lenth of the sample 
L for W = 0.2. 
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